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ABSTRACT 

We have studied the stellar population and internal structure of the core 
of the Ursa Minor dwarf spheroidal galaxy, using a sample of stars selected 
to be members based on their proper motions. In agreement with previous 
studies, we find Ursa Minor to be dominated by an old, metal-poor stellar 
population. A small number of stars with high membership probabilities lie 
redward of the red giant branch. The brightest {V ^ 18) such stars are known 
to be Carbon stars, rather than metal-rich first-ascent giants. A number of 
stars with high membership probabilities lie blueward of the red giant branch, 
and are more luminous than the horizontal branch. We speculate that these 
are post-horizontal branch stars. There may also be one or two stars in the 
post-AGB phase. Spectroscopy of the candidate post-HB and post-AGB stars 
is required to determine their nature. We recover the internal substructure in 
Ursa Minor that has been noted by several authors in the last 15 years. Using 
a variety of two- and three-dimensional statistical tests, we conclude that this 
substructure is statistically significant at the 0.005 level. There is no evidence 
that the regions of density excess have stellar populations that differ from the 
main body of Ursa Minor. The crossing time for a typical density excess is only 
~ 5 X 10^ yr. They are therefore clearly not due to intermediate age star-forming 
bursts. We conclude that they are instead due to tidal interactions between the 
Galaxy and Ursa Minor. 

Subject headings: galaxies: dwarf — galaxies: individual (Ursa Minor dSph) — 
galaxies: stellar content — galaxies: structure — Local Group 
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1. Introduction 

The discovery of the Ursa Minor dwarf spheroidal (dSph) galaxy was announced by 
Wilson (1955), although the first comments on its stellar populations were by Baade (1950). 
Ursa Minor is one of four dSph satellites of the Galaxy that were identified by direct 
examination of first epoch Palomar Sky Survey plates (Harrington & Wilson 1950; Wilson 
1955). The first color-magnitude diagram (CMD) was presented by van Agt (1967). These 
and other early photographic studies led to the conclusion that dSphs were pure population 
II systems. Their stellar populations appeared to be the same as those of Galactic globular 
clusters, but with central stellar densities lower by some six orders of magnitude. Several 
decades later, Ursa Minor remains the only dSph that appears to have a pure pop II stellar 
population (e.g., Mighell & Burke 1999; Hernandez, Gilmore & Valls-Gabaud 2000). 

Hodge (1964) performed the first structural study of Ursa Minor, and showed that it is 
a remarkably elongated system. He found e = 0.55 ± 0.10. This compares very well with 
the result of Irwin & Hatzidimitriou (1995), who found e = 0.56 ± 0.05. Hodge (1964) used 
star counts on Schmidt plates to fit the stellar distribution of Ursa Minor with an analytic 
King (1962) model. His results for the core and tidal radii are also consistent with the more 
modern values of Irwin & Hatzidimitriou (1995). 

The first paper to study the two-dimensional structure of Ursa Minor was Olszewsky 
& Aaronson (1985). This paper showed evidence for non-axisymmetric structure in the 
stellar distribution of Ursa Minor. Both the Fornax (Hodge 1961; Eskridge 1988b; Irwin & 
Hatzidimitriou 1995) and Sculptor (Eskridge 1988a; Irwin & Hatzidimitriou 1995) dSphs 
are also known to have non-axisymmetric structure in their stellar distributions. Recent 
work ( Martinez- Delgado et al. 2001) indicates that a substantial population of extra-tidal 
stars exists surrounding Ursa Minor. It is unclear if the extra-tidal population is in the form 
of tidal tails, or if it is an extended diffuse halo. The existence of this extended structure 
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argues that Ursa Minor is, in fact, not in virial equilibrium, and may have great impact on 
our understanding of the halo dwarfs. 

In this paper, we use the proper motion-selected sample of Schweitzer (1996, see also 
Schweitzer, Cudworth & Majewski 1997), discussed in §2, to study the stellar populations 
and structure of the core of Ursa Minor. In §3, we test the hypothesis that Ursa Minor is 
a pure population II system, and seek to either detect any signature of a younger stellar 
population, or to put a numerical limit on such a contribution; In §4, we apply a set of 
statistical tests to the proper motion selected sample in order to assess the significance 
of the sub-structure noted in the works above. In §5 we summarize our results, and give 
suggestions for future research. Following Mighell & Burke (1999), we adopt a distance 
modulus of (m — M)o = 19.18 ± 0.12, and foreground reddening and absorption of 
E{B -V) = 0.03 ± 0.01 and Ay = 0.09 ± 0.03, giving a linear distance of 69 ± 4 kpc. 



2. Data Acquisition and Reduction 

The data for this study are taken from Schweitzer (1996). That work provided 
membership probabilities based on proper motions for stars brighter than V ~ 20.5 (slightly 
fainter than the horizontal branch) within ~ 23' of the adopted centroid of Ursa Minor: 
0^2000 = 15^08"!8, ^2000 = +67°12'. The area surveyed is larger than the Ursa Minor core 
(rc = 15!8 ± 1(2, Irwin & Hatzidimitriou 1995). The catalog has a total of 1533 stars. Of 
these, 1060 have a non-zero membership probability. The membership- weighted sample is 
867.5 stars. 
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3. Stelleir Populations 

Figure 1 shows the V -{B-V) CMD for the full sample. The red giant branch (RGB) 
is visible from V « 17, {B - V) ^ 1.2 to V ^ 20, {B -V) ^ 0.7. The horizontal branch 
(HB) is also visible and, as is well known (e.g., van Agt 1967; Olszewsky & Aaronson 1985; 
Mighell & Burke 1999), dominated by blue HB (BHB) stars. There are also a large number 
of RR Lyrae variables in the sample. However, the CMD is contaminated by a large number 
of foreground stars. We can use our membership probabilities to produce a much cleaner 
CMD, shown in Figure 2. This figure compares well with Fig. 5 of Cudworth, Olszewsky 
& Schommer (1986), which was based on an earlier proper-motion solution with a shorter 
time baseline. In Fig. 2, the size of the symbol corresponds to the membership probability 
(P), with the smallest dots having P < 0.1, and the largest circles having P > 0.9. We 
note that the RGB is now a very well-defined feature, and that there is evidence for an 
asymptotic giant branch (AGB) as well. There are a number of stars with high membership 
probabilities in unusual places of the CMD. 

In region I, there are three stars with P > 0.5 that are substantially brighter and bluer 
than the RGB tip. The joint probability that all three of these stars are non-members 
is only w0.007. We speculate that one or more of them may be evolving across the H-R 
diagram toward the Planetary Nebula stage. None of them appear to have been studied 
spectroscopically. This seems like a project well worth undertaking. 

In region II, there are a dozen stars with P > 0.5. Although one of these stars (shown 
with a cross in Fig. 2a) is known to be a non-member (Hargreaves et al. 1994), the joint 
probability that none of them are members is < 10~^^. The members in this region would 
be post-HB or AGB stars. 

Finally, in region III, there are 17 stars with P > 0.5. Eight of these are velocity 
members based on the spectroscopy of Hargreaves et al. (1994), Armandroff, Olszweski 
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& Pryor (1995), and Olszewski, Aaronson & Hill (1995). Two are velocity non-members 
according to Shetrone, Cote & Stetson (2001b). The other seven have not been observed 
spectroscopically. The joint probability that none of these seven stars are members is 
P 5 X 10~^. Member stars in this region of the CMD are candidate Carbon stars. Indeed, 
of the eight confirmed members, five are known to be Carbon stars (Armandroff et al. 1995; 
Shetrone et al. 2001b). One is a Red Giant (Shetrone et al. 2001b), very close to, but 
slightly redward of the RGB tip. In Fig. 2a, the non-members are shown with crosses, the 
Carbon stars with squares, and the red RGB tip star with a triangle. 

In Table 1, we give the positions, magnitudes, and colors of the stars in all three regions 
of the CMD as an aid for future spectroscopic observations. We note those stars that 
are known to be velocity members based on the spectroscopy of Hargreaves et al. (1994), 
Armandroff et al. (1995), Olszewski et al. (1995), and Shetrone et al. (2001b). 

In Fig. 2b we show the ridge lines for the RGBs of a set of globular clusters, overlayed 
on the Ursa Minor CMD. We use the globular clusters M68 (Walker 1994), M55 (Lee 
1977), NGC 6752 (Cannon & Stobie 1973), NGC 362 (Harris 1982), and 47 Tuc (Hesser et 
al. 1987), with minor corrections as given in Da Costa et al. (1996), and shifted to account 
for the distance and reddening of Ursa Minor. In Figure 3, we show the probability- weighted 
CMD from Fig. 2, along with sets of theoretical isochrones from Girardi et al. (2000). As 
with the giant branch ridge lines, we have shifted the isochrones to account for the distance 
and reddening of Ursa Minor. We show a set of 14 Gyr isochrones for a range of abundances 
in Fig. 3a, and a set of isochrones with [Fe/H] = —1.7 for a range of ages in Fig. 3b. As 
shown by extensive previous work (Shetrone, Cote & Sargent 2001a, and references therein), 
Ursa Minor is dominated by an old, metal-poor stellar population. As is clear from Figs. 2b 
and 3a, there are proper-motion members of Ursa Minor that fall on the RGB loci of old, 
comparatively metal-rich stellar populations {[Fe/H] w —0.7). Recent spectroscopy of 
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some of the brighter stars in this region of the CMD shows them to be either non-members 
or Carbon Stars (Shetrone et al. 2001b). At the moment, there is no observational support 
for the existence of an old, but metal-rich stellar population in Ursa Minor. 

Figure 4 shows the B- and y-band luminosity functions (LFs) of Ursa Minor, 
incorporating our membership probabilities. The error bands are simply root-N errors of 
the unweighted number of stars. That is, any star with a non-zero membership probability 
adds one to the sum for error calculations, but adds only Pj to the LP. We apphed a Sobel 
filter (Madore & Freedman 1995; Lee, Freedman & Madore 1993) in order to search for 
structure in the LFs. We show an example of the binned differential l^-band LF, and the 
resulting Sobel-filter output in Figure 5. The RGB tip is at F «i 16.4 (or My ^ -2.9). 
There is an apparent gap at F 17.5, and the Red HB stars cause a strong feature at about 
V « 19. No other features are obvious in the Sobel-filter output. The gap at V ~ 17.5 
is apparent in the CMD. But there is no reason to expect such a gap at My ~ —1.9 in 
an old, metal-poor RGB. To investigate this further, we constructed a pure RGB F-band 
LF. We did this by selecting stars in the un-numbered region of Fig. 2a, down to V" = 19. 
In Figure 6, we show the logarithmic integrated \^-band RGB LF (solid line), along with 
Poisson error bars, computed as above (dotted hues), and a generic RGB LP with a slope 
of 0.6 (dashed hne). This shows that the apparent gap at F 17.5 is consistent with being 
a statistical fluctuation at about the 1-cr level. 

4. Tests for Substructure 
4.1. Background 

Evidence for non-axisymmetric structure in Ursa Minor was first noted by Olszewsky 
& Aaronson (1985). However, this was based on an analysis of the distribution of Ri3000 
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stars in a small (f« 3' x 5') field. Subsequent studies by Irwin & Hatzidimitriou (1995) and 
Kleyna et al. (1998) reinforce the notion that there is non-axisymmetric clumping in the 
distribution of stars in Ursa Minor, but none of these three studies has been able to make a 
firm statement on the statistical significance of these structures. 

Demers et al. (1995) rediscovered one of Olszewsky & Aaronson's stellar clumps with 
their larger-format CCD data, but did not realize this at the time, for an accumulation 
of reasons given in Battinelli & Demers (1999). BattineUi & Demers (1999) use deep, but 
very small format, HST/WFPC2 data to argue that the density excesses indicated by 
earlier work are actually due to a ring with a central depression. This is a truly remarkable 
suggestion, and warrants further study by any other means possible. 

Our data are not as deep as the Olszewsky & Aaronson (1985) CCD data, and cover 
a smaller field than the Irwin & Hatzidimitriou (1995) Schmidt plate scans. However, 
our data combine a moderate field of view (roughly 40' x 40'), with a moderate depth 
(^im ~ 20.5), and with proper-motion based membership probabilities. Thus the thorny 
problem of foreground subtraction is much less of an issue for our data than any of the 
previous data sets used for structural studies. 

In Figure 7a, we show a contour plot of the probability-weighted stellar density of our 
sample. Figure 7b shows a grey scale map of the density of non-members. We applied a 
spatial median filter of 9' to the probability-weighted membership image. The filter box-size 
is roughly the minor-axis core radius of Ursa Minor as reported by Irwin & Hatzidimitriou 
(1995). We experimented with a range of filter kernals, and found the 9' to be the best at 
removing the large-scale structure of Ursa Minor, without filtering out the regions we are 
interested in. After generating the filtered image, we then subtracted it from the raw image, 
and show the result in Figure 7c. Density excesses on either side of the centroid are now 
obvious in the median-filter subtracted data. We give the positions of these density excesses 
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in Table 2. The essential question that no previous study has resolved is the following: Are 
these substructures statistically significant? 

4.2. Analysis 

In many respects, this problem is similar to that of searching for substructure 
in clusters of galaxies. The significant distinction is that one has position and radial 
velocity information for galaxy cluster studies, whilst we have position and proper motion 
information. Our problem can be made dimensionally equivalent to that of galaxy cluster 
studies by taking membership probability as our third measurement. 

Pinkney et al. (1996) compile a set of one- two- and three-dimensional tests for 
substructure of sparse data that they have applied to galaxy cluster studies. We have 
applied a subset of these tests to our sample, using software kindly provided by Jason 
Pinkney. Specifically, we use the /3, AST, and Lee 2D tests from the two-dimensional 
substructure tests, and the Lee 3D, A, a, ftuar, and e tests from the three-dimensional 
substructure tests. We do not use the ID tests, as they search only for substructure in the 
third variable. Nor do we use the Fourier Elongation Test, as we already know that Ursa 
Minor is elongated. 

In Table 3, we show the results of the 2-D and 3-D tests. The analysis was done twice; 
once using the full data set (A^ = 1060), and once using a sub-sample of all objects with P 
within 3(7 oi P [N = 890). In practice, the 170 stars removed by the clipping procedure 
have P < 0.63. All three 2-D tests return highly significant results for substructure using 
both the clipped and unchpped samples. Amongst the 3-D tests, the Lee 3-D test returns 
highly significant results for substructure using both the clipped and undipped samples. 
The A, a, and a^ar tests consistently give no positive results for substructure, and the e 
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test argues for substructure only when applied to the clipped data set. 

The j3 test (West, Oemler & Dekel 1988) returns evidence for statistically significant 
subclustering in Ursa Minor on both the clipped and undipped samples. The (5 statistic 
is a measure of deviation from mirror symmetry. Pinkney et al. (1996) note that the (5 
statistic does not measure deviation from circular symmetry, so our result is not simply due 
to the known elongation of Ursa Minor. Rather, this test argues that the distribution of 
stars about the galaxy center is clumpy, with the clumps not being distributed according to 
mirror symmetry. 

The angular seperation test, or AST (West et al. 1988) also returns evidence for 
statistically significant subclustering in Ursa Minor on both the clipped and undipped 
samples. This test measures the excess of small-angle seperations above that expected 
due to Poisson noise. West et al. (1988) examined the behavior of the AST on elongated 
structures, and found that false positives did not occur for ellipticities up to e ~ 0.5 As 
Ursa Minor has e = 0.56 ± 0.05 (Irwin & Hatzidimitriou 1995), this is unlikely to be the 
cause of the positive signal we see for the AST. 

The Lee 2-D test (Lee 1979) also returns evidence for statistically significant 
subclustering in Ursa Minor on both the clipped and undipped samples. This test measures 
the significance of seperating the data set into two subclusters. As with the (3 test, an 
elongated but smooth distribution will not cause the Lee statistic to report a false positive. 

The Lee 3-D test (Fitchett & Webster 1987) is an extension of the Lee 2-D test that 
incorporates variation in the third dimension. We note that if the spatial subclusters have 
essentially the same distribution in the third variable (membership probability, in the case 
of our Ursa Minor data), the Lee 3-D test should give the same result as the Lee 2-D test. 
This is the result we find for our Ursa Minor data. 
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Of the other 3-D tests, the A, a, and a^ar tests consistently give no positive results for 
substructure. The A statistic (Dressier & Shectman 1988) looks for significant variations in 
the third dimension for spatial clumps. As our third dimension is membership probability, 
the null result we obtain simply says that the spatial subclusters (if any) are part of 
Ursa Minor. That is, the stars in the spatial subclusters have the same distribution of 
membership probability as do the rest of the Ursa Minor stars. The a (West & Bothun 
1990), and a^ar (Pinkney et al. 1996) tests both measure the dispersions of the third 
dimension for spatial clusters, and thus, like the A statistic, would not be expected to show 
significant subclustering in the present case. 

This brings us to the e test (Bird 1993). In this context, e is a measure of the 
average number of stars per nearest neighbor group. Paradoxically, this test gives strong 
evidence for subclustering on the clipped data set, but no evidence for subsclustering on the 
undipped data set. This appears to be due to the diluting effects of the low membership 
probability stars on the test statistic (see discussion in Pinkney et al. 1996). When only the 
high membership probability stars are retained, the dispersion in membership probability 
for any spatial subcluster will be much smaller, leading to a more significant result for the e 
test. 

The essential result of this set of tests is that we confirm the previous evidence for 
subtructure in Ursa Minor and can now quantify the existence of this substructure with a 
statistical confidence of > 0.995 for all well-posed tests of that structure. The distribution 
of member stars in Ursa Minor is not well described by a smooth distribution. In §5, we 
discuss the imphcations of this for the evolution of Ursa Minor. 
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5. Summciry and Discussion 

We have analysed a proper-motion selected sample of stars in the Ursa Minor dwarf 
spheroidal galaxy in order to study the stellar populations, and structure of this object. 
The data are taken from Schweitzer (1996), and include stars within ~23' of the Ursa Minor 
centroid, down to a limiting magnitude of y ~ 20.5. The membership-weighted sample is 
867.5 stars from a total sample of 1533. 

The CMD for our sample compares very well with previous bright-star CMDs for 
Ursa Minor (van Agt 1967; Cudworth et al. 1986), showing the dominant old, metal-poor 
stellar population. Comparison with globular cluster RGB ridge lines and model isochrones 
indicates a metallicity of [Fe/ H] w —1.8 to —2, consistent with the spectroscopic result of 
[Fe/H\ — — 1.90±0.11 from Shetrone et al. (2001a). There are a number of high- probability 
members in unusual parts of the CMD. We provide position and magnitude data for these 
stars as an aid for future spectroscopic observations. The B- and l^-band LFs of our sample 
are consistent with a single, dominant, old, metal-poor stellar population. 

We have applied a number of tests for substructure to our sample in order to evaluate 
the statistical significance of the well-known clumps of stars in Ursa Minor (Olszewsky 
& Aaronson 1985; Irwin & Hatzidimitriou 1995; Demers et al. 1995; Kleyna et al. 1998; 
Battinelli & Demers 1999). These tests are drawn from the study of Pinkney et al. (1996), 
who compiled a set of one- two- and three-dimensional tests for substructure in sparse data 
for galaxy cluster studies. When the detailed assumptions and behaviour of these tests 
are taken into account, we conclude that the hyphothesis that Ursa Minor has statistically 
significant substructure is strongly supported (see Table 3) . We also employed a test devised 
by Mighell (2001) to search for structure in sparse 2-D Poisson-distributed data. The test 
is a modification of the statistic. Given a model, the test reports the significance of the 
deviation of the data from the model. As the test is well-defined for integers, we have used 
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the sample of all stars with P > 0.9 of being members of Ursa Minor. Our model is a simple 
elliptical profile fit to these data, with ellipticity, and position angle as free parameters. 
The test rejects the null hypothesis that the data follow the model at the 95% confidence 
level. Thus, this test is consistent with the results from the clustering statistics discussed in 
§4, and indicates that the distribution of stars in Ursa Minor is not well-fit with a smooth 
model. 

We do not see evidence of the structure reported by BattineUi & Demers (1999). The 
lack of agreement between our study and theirs is likely due to the very different strengths 
and weaknesses of the data for the two studies. The BattineUi & Demers (1999) study 
is based on very deep data covering a very small field of view (a single WFPC2 frame), 
whereas our data are very shallow, but cover a much larger field of view. A proper test of 
the results of BattineUi & Demers (1999) would require much deeper wide-field data than 
we have at our disposal. 

The angular size of the lumps is roughly 3' (see Figure 7c). This is about 20% of 
the size of the core radius. The measured velocity dispersion of Ursa Minor is roughly 10 
km/sec (10.4 ± 0.9 km/sec according to Armandroff et al. 1995). For our adopted distance 
of 69 ± 4 kpc, this implies a lump crossing time of only ~ 5 x 10^ yr. However, neither the 
CMD, nor the \^-band LF of stars in the lumps shows any evidence for differences in the 
stellar populations of the lumps and the bulk of Ursa Minor. The lumps are clearly not 
due to young (or even intermediate-age) star formation events. We speculate that they are 
instead due to tidal stretching from the gravitational interaction between the Galaxy and 
Ursa Minor. Features like the Ursa Minor lumps appear in some numerical simulations of 
such interactions (e.g., Kroupa 1997; Klessen & Kroupa 1998). However, the evidence for 
such dynamically induced structure and for extra-tidal stars associated with Ursa Minor 
( Martinez- Delgado et al. 2001) is unlikely to provide an escape for the large mass to light 
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ratio for Ursa Minor that is implied by its velocity dispersion, as numerical studies indicate 
that the observed velocity dispersions of disrupting dwarfs do not significantly exceed their 
viral values until just before their complete dispersal (e.g., Oh, Lin & Aarseth 1995; Piatek 
& Pryor 1995). 

It would be very interesting to obtain spectra of the remaining high-probability 
member stars in the outlying regions of the CMD. The stars in region I are bright enough 
{V < 16.2) that spectroscopy with a 4m class telescope should be able to determine their 
membership, and their nature. All the high-probability stars in regions 11 and III that have 
not yet been observed spectroscopically are fainter than V ^ 18, and thus become difficult 
targets for 4m-class spectroscopy. Determining the membership and nature of these stars 
would be an excellent project for the HET, the upgraded MMT, or the Keck telescope, 
or, in the future, for the Gemini telescope. Very deep photometry of Ursa Minor, both 
in the lumps and in the general field, would provide a much better probe of the stellar 
populations than the data currently available. The main problem with the WFPC2 data 
is that the field of view is so small. Observations of a number of fields in the halo dwarf 
spheroidals with the MOSAIC Imager on the Mayall and Blanco 4m telescopes, and with 
the forethcoming Advanced Camera for Surveys will put our understanding of the details 
of the star formation histories of the halo dwarfs on a much better observational footing 
than is currently possible. Finally, there is now clear evidence that roughly half of the halo 
dSphs have internal substructure. The lesson of the Sagittarius dSphs is that interactions 
with the Galaxy can cause distortions in the distributions of stars in the halo dwarfs. The 
current observational situation presents a challenge to theorists that has the potential to 
lead to dramatic improvements in our understanding of the role of tidal stripping in the 
dissolution of dwarf systems, and the building of the Galactic halo. 
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Table 1 - CMD Outliers with P > 0.5 
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Notes: 

1: Numbers from Cudworth et al. (1986). 
2: Numbers from M. Irwin (unpublished). 

3: Velocity membership status from Hargreaves et al. (1994), Armandroff et al. (1995), 
Olszewski ct al. (1995), or Shetrone et al.(2001). C: Carbon stars according to Armandroff 
et al. (1995), or Shetrone et al.(2001). R: RGB star according to Shetrone et al.(2001). 
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Table 2 - Lump Positions 



Number 


RA (2000) 

hhmmss 


Dec (2000) 

o o ' ' " " 


1 


15 08 00 


67 06 30 


2 


15 08 30 


67 10 20 


3 


15 08 40 


67 12 30 


4 


15 09 10 


67 11 40 


5 


15 10 10 


67 21 50 


6 


15 10 30 


67 19 50 



Table 3 - Substructure Test Results 
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Fig. l.— V-{B-V) CMD of the full Ursa Minor sample. 

Fig. 2. — V - {B — V) CMD incorporating probability weights. The smallest sohd dots 
have membership probabilities of P < 0.1. The largest open circles have membership 
probabilities of P > 0.9. Stars with P > 0.5 are shown as open symbols. The solid lines 
demark regions of the CMD, labelled with roman numerals I-III, from which the stars listed 
in Table 1 were drawn, a). The squares show spectroscopically confirmed member Carbon 
stars, the triangle shows the spectroscopically confirmed member RGB star, the crosses show 
the spectroscopically confirmed non- members, b). The curves show globular cluster RGB 
ridge hues taken from the literature. The clusters shown are 47 Tuc (dotted line), NGC 362 
(short dashed line), NGC 6752 (long dashed line), M55 (dot - short dash hne), and M 68 
(dot - long dash line). 

Fig. 3.— V - {B -V) CMD from figure 2, with isochrones from Girardi et al. (2000) 
overlayed. a) 14 Gyr isochrones with [Fe/H] = —1.7 (solid line), -1.3 (dotted line), -0.7 
(short-dashed line), -0.4 (long-dashed hne), (short-dashed-dotted hne) and -1-0.2 (long- 
dashed-dotted line), b) [Fe/H] — —1.7 isochrones with ages of 18 Gyr (solid line), 14 Gyr 

(dotted line), 11 Gyr (short-dashed line), 9 Gyr (long-dashed line), 7 Gyr (short-dashed- 
dotted line) and 5.6 Gyr ((long-dashed-dotted line). 

Fig. 4. — Logarithmic probability-weighted integrated LF of Ursa Minor. The solid lines are 
the LFs, and the dashed hues represent root-N error bounds of the unweighted number of 
stars, a) P-band LF. b) F-band LF. 

Fig. 5. — a) Differential F-band LF with 0.2 magnitude binning, b) Sobel-filter output of 
the data shown in a). 

Fig. 6. — Logarithmic F-band integrated LF of the RGB of Ursa Minor. The sohd hne is 
the LF, the dashed lines represent root-N error bounds of the unweighted number of stars. 
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the dashed line is a model LF with a power law index of 0.6. 

Fig. 7. — Density maps of stars in the Ursa Minor sample. The origin is the adopted center of 
Ursa Minor. North is up, East to the left, a) Contour plot of the probability- weighted density 
of Ursa Minor stars binned into a 1' x 1' image, b) A grey-scale plot of the non-members 
at the same resolution, c) A contour plot of the residuals of the membership-weighted array 
minus a median-filtered version of the array, with a 9' filter box. 
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